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Abstract 

A general variational principle of classical fields with a Lagrangian 
containing field quantity and its derivatives of up to the N-order is pre- 
sented. Noether's theorem is derived. The generalized Hamilton-Jacobi's 
equation for the Hamilton's principal functional is obtained. These results 
are surprisingly in great harmony with each other. They will be applied 
to general relativity in the subsequent articles, especially the generalized 
Noether's theorem will be applied to the problem of conservation and 
non-conservation in curved spacetime. 

1 Introduction 

The aim of this series of articles is to explore the conservation and non-conservation 
in curved spacetime, especially to explore the difficulty of energy- momentum 
conservation in general relativity and the gravitational energy-momentum. We 
start with presenting a general variational principle for classical fields with a 
Lagrangian containing the field quantity and its derivatives of up to the N-th 
order (part I), Then the general results from part I are applied to general rela- 
tivity, especially the generalized Noether's theorem is applied to the problem of 
conservation and non-conservation in general relativity (part II). The last part 
(part III) is devoted to the difficulty of conservation of energy-momentum in 
curved spacetime and to the problem whether the metric field carries energy- 
momentum or not. 

The developments of modern physics, such as the founding of statistical 
mechanics and quantum mechanics, have proved that the variational princi- 
ple approach to dynamics is not only an alternative and equivalent version to 
the naive, intuitive approach, but also yields deeper insights into the underly- 
ing physics. For instance, it is hard to imagine that the statistical mechanics 
could have been established without using the concepts of phase space, and the 
quantum mechanics could have been established without using the concept of 
Hamiltonian. Therefore, we will found our argument on a general variational 



principle for classical field. It might be for the same reason, soon after Ein- 
stein proposed his general theory of relativity, Hilbcrt made the first attempt 
to get Einstein's equation by using the least action principle. The Lagrangian 
being used for vacuum Einstein's equation, (16irG)~ 1 R, is the only indepen- 
dent scalar constructed in terms of the metric field and its derivatives of no 
higher than the second order. However, because the Ricci scalar curvature R 
contains the second order derivatives of the metric field g^ u (x), which is now 
the dynamic variable, the least action priciplc for Lagrangians containing only 
the field quantity and its first order derivatives does not lead to Einstein's 
field equation. The generally accepted solution to this difficulty is adding the 
Gibbons-Hawking boundary term to the Hilbert action and keeping the least 
action principle unchanged [1]. But there is another solution to this difficulty, 
which is adopted in the present paper. The least action principle will be restated 
and the Hilbert action will still be used for the vacuum Einstein's equation. In 
order to show this is proper and natural, we will consider classical fields with a 
Lagrangian containing the field quantity and its derivatives of up to the N-th 
order. In our opinion, acting at a distance is not acceptable, so we assume that 
the Lagrangian does not contain the integral of the field quantity. In section 
2, a general Lagrangian formalism for classical fields is presented. In section 3, 
the Hamiltonian formalism is discussed. In section 4, the Noether's theorem is 
derived and the conservation law due to the " coordinate shift" invariance is es- 
tablished. The generalized Hamilton- Jacobi's equation is obtained in section 5. 
All the results obtained above are in great harmony with each other and apply 
to various classical fields, say, those with Galilean covariance, with Lorentzian 
covariance, with general covariance, or without such covariance. Part I finishes 
with a remark. In part II, this general variational principle of classical fields 
developed in part I is applied to general relativity and quite a few conserved 
quantities corresponding to the coordinate " shift" invariance, coordinate " rota- 
tion" invariance etc. are found. And the properties of these conservation laws 
are discussed. In part III, after some general consideration, the introducing of 
gravitational energy-momentum is reviewed. The difficulties of conservation of 
energy- momentum in general relativity are explored by using Noether's theorem 
and observations from geometry. It is pointed out that the metric field docs not 
carry energy-momentum, and the law of conservation of energy-momentum no 
longer holds in curved spacetimc. 



2 Lagrangian formulation of classical fields 

Suppose that our spacetime M is a smooth manifold which is differentially 
homcomorphic to K 4 . Choose a chart (M,tp), and denote by (x° , x 1 , x 2 , x 3 ) 
the corresponding coordinates. Suppose the action over any spacetime region 
£1 C M of the classical field {$ a (x)} is 

A= [ d 4 xL(x,<$>(x),d$(x),d 2 $(x),...,d N $(x)) =A[$] (I) 
Jn 
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where the Lagrangian L is a function of the spacetimc coordinate, the field 
{$ a (x)} and its derivatives of no higher than the N-ih order. L docs not 
contain the integral of {$ a (x)}, since acting at a distance is not acceptable. 
In order to develop a general variational principle for all locally interacting 
classical fields, here we suppose that L can manifestly contain the spacetime 
coordinates, and the Galilean invariance, Lorentzian invariance or the general 
invariance are not assumed for the time being. For the sake of simplicity, we 
have assumed that our spacetime manifold is (3 + l)-dimcnsional. However, our 
presentation has nothing to do with the spacetime dimensionality. It still holds 
for an (n + l)-dimcnsional spacetime. 

Consider the difference between actions over f2 of two possible movements 
close to each other. Using integration by parts and Stokes theorem, one gets 



6 A[^] = f d 4 x6^ a (x)[-^- - d Xl + _ 

J n d<$> a {x) 'ddx^aix) 



dL 



+ (-l) N d Xl ---d XN - T ]+ / ds x [B^53> a {x)+B^5d v ^ a {x) 

dd Xl ---d XN <P a (x) J on 

+ B aXvil/2 5d Vl d V2 $a(x) + ■■■ + B aXvi "" N - 1 5d Ul ■ ■ ■ d, N _^ a (x)}, (2) 
where the Greek indices go through 0,1,2,3, and 



r>a\ _ d L fl dL dL 

dd x $ a (x) ^ddxd^ a {x)^ ^ ^dd x d^d^ a {x) 



dd x d^---d^ N _^ a { x y 

B ^ = d ± d — + d d — + ■ 

ddxd vi * a (x) 111 ddxd^dMx) ^ ^ 112 ddxd^dMx) 

+ (-l) N - 2 d — 



dd x d Vl ■ ■ ■ dv N _ 2 $a{x) Ml dd x d^d Ul ■ ■ ■ d UN _ 2 $ a (x) ' 
dT 

dd x d vl ---d VN _^ a { x y y> 

Equation (2) suggests the least action principle reads as follows. 
For any spacetime region f2, among all possible movements in f2 with the 
same boundary condition 

St>\ d n = 0, 5d$\ dn = 0, . . . , 6d N - 1 <P\ dn = 0, (4) 

the real movement corresponds to the stationary value of the action over ft. 
Combining eqns.(2), (4), one obtains the field equation ( Euler-Lagrange 
equation ) satisfied by the real movement 
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9 A 1 TTT r -^ 7 - T +- •••+(-!) 5 Ai ---5 An — — - = 0. 



<y* (x) 9$ a (x) A1 dd Al $ a (x)' v ' A1 AN a9 Al ---a Ajv $ a ( x ) 

(5) 

3 Hamiltonian formulation of classical fields 

To formulate the Hamiltonian formalism of classic fields, one needs to specify a 
reference coordinate system such that the hyper-surfaces, S^o,of 

constant £° are spacelike Cauchy hyper-surfaces and the curves, t-g ,of constant 

£ are timelike world lines of the observer at £ . We will consider the state of 
the field on hyper-surfaces, X^o, and investigate the change of state (evolution) 
with £°. We will observe the state on S^o and the evolution with £° from any 
reference coordinate system in the same way. As has been pointed out [2], the 
properly formulated Hamiltonian formalism is compatible with all dynamic sys- 
tems, Galilean invariant, Lorentzian invariant, general invariant and so on. The 
invariance is the heritage from the Lagrangian being used. The 3 + 1 decom- 
position of spacetime proposed above is more general than the one generally 
accepted in General relativity. The latter relies on the unknown dynamical 
variable, the metric field. Suppose the Lagrangian functional is 

A = J d 3 £L& 0*(O. • ■ • - d N m) = A[£°, $| £0 , d $\ e , . . . , d^\ e ], 

(6) 

Consider the difference between Lagrangian functional of two states close to 
each other. Using the general formula (2), which is independent of the dimen- 
sionality, one gets 



fir 



■ r 9L dL 

BT 

+ (-i)»-c8A-a.- a8ti ^ a j 3 g iw.(fl 
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dL 



+ [Kf5do<S> a (0 + KT kl 5d kl d^ a {0 + Kf k ^5d kl d k2 do$ a (0 

+ --- + K a 1 lkl - kN - 2 8d kl ---d kN _ 2 d $ a (0} 

+ [K^5d^ a (0 + Kr kl 5d kl d^ a (0 + Kr klk2 Sd kl d k2 d^ a (0 

+ ... + K? kl - k »- 3 6d kl ■ ■■d kN _ 3 d 2 <f a (0] 

+ ■■■ + [K$- 2 5d?- 2 * a (Z) + K a N %Sd k X^a(0] + Kti-iSd*- 1 ^)}, 



(7) 

where the domain of integration is M 3 , the latin indices go through 1,2, 3, and 



dL dL 8L 

ddi* a (Q 31 ddidMt) + 31 32 ddid^Mt) 



A-aifel = d ± Q. Q ± +Q.Q. Q ± + . . 



+ (-i) N - 2 d n ---d 0N _ 2 



dL 



dd l d n ■■■d lN _ 2 d kl ^ a (0' 



K aik^-k N - 2 ^ q dL 



j^aiki---kw -l 



dd t d kl ■ ■■d ktl _Mt) 31 dd l d n d kl ■ ■ -d kN _MO ' 

dL 

dd t d kl ---d kN _MO' 



K^-c 1 91 c'd dL c '3 d dL 



■ddidoQait) 3 31 ddidfrdoQaiO 4 31 32 dd t d 3l d n do* a (Z) 

+ " ' + ( - 1)JV " 2cjv ^ ' ' ' nd.n, ...«, x A ,<KW 

x « fel = l ^ _ 1, d£ , -lo o , . 

+ (-l) N - 3 c],d n ---d jN _ 3 " L 



ddid n ■ ■ ■ d jN _ 3 d kl d $ a (£) 
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,ajfci---fc]v_3 _ 1 

11 Ar - 1 9^ • • • d kN _ a do* a {0 N 31 dd t d n d kl ■ ■ ■ d kN _ 3 d Q $ a (t) ' 



T^aik 1 ---k N - 2 _ A dL 



dd l d kl ■ ■ -dkn^doQait)' 



K ai _ c 2 9L _ dL dL 

3 dd^ a (0 c ^dd,d n d^ a (0 + 5 31 32 dd t d n d J2 d^ a (0 

R-a«si = r 2 dL 2a. dL , 2a. a. ^ , . 

, C_1\JV-4 2 o o 

1 J ^^%...^_ 4 a fcl 9 2 <i> Q (e)"-" 

j^aifei - fcjv-4 _ 2 2 a 

-"■2 -^-1"^"^ a xZxrTt\ C N°h 



1 ddid kl ■ ■ ■ d kN _ 4 d^ a (0 N 31 ddid h d kl • • • 9 fejv _ 4 9 2 $ a (C) ' 



i\ 2 — c N 



dd t d kl ---d kN _ 3 di$> a (zy-- 



K ai _ N-2 dL r N - 2 8 9L 



Hence 



ddidS'-^aiO ' N 31 ddidjX' 2 **®' 
K ai kl _ N-2 ^ K ai _ „N -\ 

dd l d kl d $ a (£) 9(9,9™ i $ a (^) 



(5A „ <9L 



5*»(0 a*„(0 n ddiMO 

FIT 



SA dL dL 



+ ---- + (-i) N - 2 c%d il ..-d i , i)L 



dd il ---d iN _ 2 dis> a (ay 
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SA dL 



(9) 



Consider the difference between actions of two possible movements close to 
each other. 

6A= /V°A[^o,d <%o,...,^<%o]= [ d? [ d 3 CL($(0,dm,---,d N m) 

J to J t^ J 

+ ■■■+ sdN* a ® S9 °* am + I] ^ J ^U*™'® + K a ^5d kl <P a (0 

+ K aik ^Sd kl dMO + ■■■ + K aikl " k »-*6d kl ■ ■ ■ d kN _MO + K aik i- k »5d kl ■ ■ -0 fcw *„(£)] 

+ [K?5do$ a (Q + K a ^6d kl dM0 + Kf k ^5d kl d k2 d $ a (0 + ■ ■ ■ 

+ Kf kl - kN -'5d kl ■ ■ ■ d kN _ 2 d $ a (0] + [K?5%$ a (0 + K^S8 kl d^ a (0 

+ K a 2 ^8d kl d k2 dl<s> a {0 + ■■■ + Kf kl - kN ~ 3 6d kl ■ ■ ■ d kN _ 3 d 2 M0] + ■■■ 

+ [*#-2^ _2 *«(0 + K a ^%5d k X^a(0] + K^Sd?- 1 ^)} 



+ B aX ^Sd Ul d u MO + ■■■ + B aXvi "-"'- 1 6d Vl ■ ■ ■ d UN _^ a (0 
+ B aXv *"'"'6d Vl ---d v MQ] (10) 
Using the least action principle, one re-obtains the Euler- Lagrange equation 

^ A o ^ A | n2 SA N N 5 A 

Noting eqn.(9), one easily sees that eqns.(ll) and (5) are exactly the same. 
Let 



(12) 
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One easily gets 



SH = J rf 3 e[9o$a(?)K(0+5 2 <fa(0^2(?) + - ' ■+^*a(0<^(0-£^*MO] 

-J dai{[K ai 5<S> a (® +K a ^6dMO + K a ^ k *Sd kl d k2 <t> a (0 

+ --- + K aik i- k »-i5d kl ---d kN _ 1 $ a (0] 

+ [KfSdo<S> a (Z) + Kf^6d kl d ^ a (O + K« iklk2 sd kl d k2 d a <s> a (0 

+ ... + K? kl - kN - 2 5d kl ■ ■ ■ d kN _ 2 doMO] 

+ [K?6a$* a (0 + Kr kl sd kl d 2 MO + Kr klk2 Sd kl d k2 d 2 MO 

+ ... + K« lkl - kN - 3 sd kl ■ ■ ■ d kN _ 3 d 2 $ a (0} 

+ ■■■ + [K#_ 2 5d»- 2 * a (Q + K%%5d k X- 2 $a((;)] 
+ K%_ 1 5d^- 1 <i> a (0}. (13) 
This suggests that the Hamiltonian H is a functional of {<& , n", . . . , n^} 

H = H[?,*,K U n 2 ,...,n N ], (14) 
5H SA SH o i rc , SH SH N 

s^) = -M^y^m= s^m = d ^ ■ ■ ■ ' s^T) = d ° ^ 



(15) 



and 



A - L + ^ S(S) + ' ' ' + <16) 

From the Eulcr-Lagrange equation (11), one gets 

SH 

fane d 2 n a 2 (0 - + ■■■- (-i) N d^%(0 = -j^ y (17) 

Eqn.(17) and eqn.(18) 

*••<« = J^)' a »*««> " a "* ae) = 4®' (18 ' 

constitute the canonical equations. Note that when N = 1 ( all pre-G.R. 
field theories belong to this case ), canonical equations (17), (18) read 
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And one has 

±H = d«H. (20) 
When N = 2 ( G.R. is this case ), the canonical equations read 

SH 



9o<(0-«(0 



And one has 



^{H - J d 3 £[d $ a (£)d ^(£)]} - d H. (22) 

4 Noether's theorem 

4.1 Proof of Noether's theorem for Lagrangians contain- 
ing up to iV-th derivatives of field 

Nowwehaveto deal with two kinds of derivatives of L(x, 3>(x), d$(x), <9 2 $(x), . . . ,d N $(x)) 
with respect to coordinates, d„L = dLjdx" and (5 a L = dL/(5x a , relating to each 
other through the following equation, 

9L 9L dL . dL 

ox' 7 ox a oq> a {x) oox 1 9 a (x) 

+ ■•■+»» » *> ( , 9 a d Xl ---dxMx) (23) 

od\ l ■ ■■d\ N $a(x) 

Theorem 1 If the action of classical fields over every spacetime region fi re- 
mains unchanged under the following r— parameter family of infinitesimal trans- 
formation of coordinates and fields 

x x — > x x = x x + 5x x , 

*a(s) $a(x) = * a {x) + 5$ a { X ), (24) 

then there exist r conserved quantities. 
Proof. From eqn. (24) one has 
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Sd 4 x = (d a Sx a )d 4 x 1 
Sd x = -(dx8x°)d a , 
5[d Xl <f> a (x)} = (Sd Xl )^ a {x) + d Xl 5$ a (x) = d Xl 6$ a {x) - d a $ a (x)d Xl 6x°, 
6[d Xl d X2 <£>a(x)} = d Xl d X2 6$ a (x) - d X2 d a $ a (x)d Xl 5x a - d Xl d a ^ a (x)d X2 Sx a 

-d„$ a (x)d Xl d X2 5x°, 



5{d Xl d X2 ---d XN <P a (x)} = d Xl d X2 ---d XN 5$ a (x)- E d a d Xl ■■■a A( ■ ■ ■ d x Mx)d Xt 8x ,T 

l<i<N 

- E d a d Xl ---d Xi ---d Xj --- d XN $ a (x)d Xi d Xj 8x° 

l<i<,j<N 

E 9 <T d Xl ■■■d Xi ---d Xj ---d Xk --- d XN $ a (x)d Xt d X] d Xk 5x a - {!•) 

l<i<,j<k<N 

-J2 d <r d ^a(x)d Xl ■ ■■d Xk ■ ■■d XN Sx a 
k — 

-d^ a (x)d Xl ■■■d XN 5x° (25) 
Substutute eqns.(24) and (25) into the following equation 

r r ar fir f)j 

+ -" + ^ — --^••■a A „$ (1 (i)] 

99 Al •••a Alv $ a (x) 



one gets 

dL 

</*«(*) WA1 aa Al * a (x) ' ' v ^ y 9a Al ---a Ajv $ a ( 2; ) J 



7o L (9$ a (a;) 1 99 Al $ a (a;) 
(5$„(a:) - 9 ff $ a (z)<5z CT ) + / d 4 zd A [I^&r ff + B aX (61> a (x) - d^ a (x)5x a ) 
+ B aXl ^d Vl {5<5> a {x) - d a $> a {x)8x°) + B aXv ^d„,d V2 {5$ a {x) - d a <5> a {x)5x°) 

+ ... + Ba^i-tN-i^ . ..dy^SQaix) - d^ a (x)5x")} (26) 

The first integral at rhs vanishes for real movement, hence the second integral 
does too. One gets the following equation due to the arbitrariness of tt. 

= d x [LS x 6x a + B aX (S$ a (x) - d^ a (x)8x a ) + B aXvi d Vl (6$ a (x) - d„$ a (x)6x°) 
+ B aXv ^d Vl d V2 (5$ a (x) - d a <S> a (x)5x°) + ■■■ 

+ B aXvi '" VN ~ 1 d Vl ■ ■ ■ d^^aix) - d a <5> a {x)5x°)] (27) 

Noting that both 5x a and 5<& a (x) depend on r real parameters, one can consider 
eqn.(27) as r conservation laws. ■ 
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4.2 "Conservation law due to "coordinate shift" invari- 



ance 

In this subsections, we restrict our discussion to Lagrangians which do not 
manifestly contain coordinates and is invariant under "coordinate shift". In 
this case, the action (1) remains unchanged under the following "coordinate 
shift" . 

Sx a = e a , 8<S> a {x) = 0. (28) 

In this case, eqn.(27) reads 

dxr x a = 0, (29) 

where 



r x = B aX d a $ a {x) + B aXvi d Vl d c7 <& a {x) + B aXuilJ2 d Vl d V2 d a <& a {x) 

+ ■■■ + B aXvi "' VN - 1 d Vl ■ ■ ■ d^^OMx) - LS X (30) 

is usually called energy- momentum tensor. Notice that "coordinate shift" cqn.(28) 
is not an invariant concept under general coordinate transformation. This is eas- 
ily seen from the active viewpoint of transformation. This explains why in 
not a tensor under general coordinate transformation. We will get back to this 
problem later. 

5 Hamilton's principal functional and Hamilton- 
Jacobi's equation 

Let us consider the difference between actions over spacetime region of two 
real movements close to each other. Using eqns.(2) and (5), one gets, for real 
movements 



6A[$] = / ds x [B aX 8$ a (x) + B aXvi 8d Ul $ a (x) + B aXv ^5d Vl d U2 $ a {x) 
Jen 

+ ■■■ + B aX ^- v »-^Sd Vl ■ ■ ■ d vN -^ a {x)] (31) 

From eqn.(31), one sees that the action over a spacetime region Q, of a real 
movement is determined by the closed hyper-surface dtt, and <&|an, d<& | gn, ... , 
d N ~ 1 $>\gn. It will be called the generalized Hamilton's principal functional and 
denoted by 

S = S[dQ, $| an , 9$|ao, • • • , ^"^lan] (32) 

Re- write eqn.(31) as 
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SS= f ds x [B aX S^ a (x) + B aX ^Sd^ a (x) + B aXv ^5d Vl d V2 <S> a {x) 
+ ■■■ + B aX ^- v »-^5d Vl ■ ■■d vN - 1 * a {x)\ (33) 

Note that when <j> a |dfi is given, only one of the four derivatives d\<& a \dn (A = 
0, 1, 2, 3) is independent; when d\<& a \dn is given, only one of the four derivatives 
d^d\$ a \gn (fi = 0, 1, 2, 3) is independent; and so on. Thus for a given suffix a, 
only AT items from $ a |da,<9Ai$o|dn, ■■■,d\ 1 - ■■d\ N _ 1 $ a \dn (Aj = 0, 1,2,3) are 
independent. 

In order to formulate the generalized Hamilton- Jacobi's equation, one needs 
a new type of functional derivative. 

Definition 2 Let E be a hypersurface in spacetime M , a function defined on 
M , and F — F[T,, a functional of E and ^>\s. The functional derivatives 
are defined as follows. If the variation of F can be written as 

5F[E, *| E ] = J ds x {y[S, *| E , x) x 6V(x) + Z[E, *| E , x) x 5HF(x)} (34) 

then Y[S,^\^,x) x and Z[E, vE^s, x) x are called the functional derivative of F 
with respect to E M (x) and ^(x), and denoted by 

respectively. 

Hence we have 



SF ^ = L ds > f (^) ) A 6mx) + (^M ) A 



(36) 



The hypersurface E is given by the parameter equation 

x^ 1 = E"(0\0 2 ,0 3 ) (37) 

The (5E^(x) in cqn.(33) is 

<ys"(x) = w{e\e 2 ,e 3 ) - w(e\e 2 ,e 3 ). (38) 

Now, from eqn.(33) we have 



SS \ A „„i / (5S* V „„i„. / SS 



5$ a (x)J ~ '\Sd^ a (x)J ~ ' ''"'\5d Vl ---d VN _ 1 * a (x)J 

(39) 

Follow the evolution of one real movement and observe the change of its action. 
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SS= I ds x L8 x JY, a {x) 

From eqns.(30), (39) and (40), we get the generalized Hamilton- Jacobi's equa- 
tion. 

(«^) A+ ^ = »- < 41 > 

Remark 3 So far we have presented a general variational priciple for clas- 
sical fields. The only postulate made in this formalism is the least action 
principle. This formalism applies to all the classical fields {<fr a (x)} with a 
Lagrangian L(x, &(x), <9$(x), d 2 &(x), . . . , d N $(x)), say, Newtonian fluid me- 
chanics, Maxwell's electromagnetic field, general relativity, etc. The specific 
symmtries and covariance of a classical field are the heritage from the La- 
grangian, not from this general formalism. This formalism yields manifestly 
Galilean (Lorentzian, general) covariant field theory when the inputted Lagrangian 
is Galilean (Lorentzian, general) covariant. It is worth noting that all the re- 
sults obtained above, are in great harmony with each other. We will apply this 
general variational priciple to general relativity, especially apply the general- 
ized Noether's theorem to the long standing problem, conservation and non- 
conservation in curved spacetime in part II and part III. 

Acknowledgement 4 I am grateful to Prof. Zhanyue Zhao, Prof. Shihao 
Chen and Prof. Xiaoning Wu for helpful discussions. 
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